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1. Introduction 

The discussion on the stability of quasi-crystals has not resulted in a general consensus 
yet. Even very recently, arguments against |I| and in favour of the random tiling 
scenario have appeared in the literature. It has also been suggested that the entropically 
stabilised state results from quasi-crystal growth From the point of view of statistical 
mechanics random tiling models are very interesting, not in the least because some of 
them allow for an exact solution. This means that one is able to derive exact expressions 
for the entropy and other thermodynamic quantities. In this paper we present a random 
tiling with a decagonal phase which is solvable by the Bethe Ansatz method, very 
much in analogy with the dodecagonal square-triangle and the octagonal rectangle- 
triangle tilings 0, 0, 0. The random tiling with a decagonal phase is the ensemble 
of tilings of the plane by rectangles and isosceles triangles with sides of length 1 and 
/ = 2sin(7r/5) = a/2 + r/r, where r = {y/b + l)/2 is the golden mean. Cockayne 
devised an inflation rule for this set of tiles that constructs a tiling with decagonal 
symmetry. This tiling corresponds to a maximally dense decagonal disc packing under 
the condition that nearest neighbour vectors are limited to certain directions. The 
random tiling with additional constraints has been studied by Oxborrow and Mihalkovic 
to model d-AlPdMn. More recently. Roth and Henley |T0| used the rectangle-triangle 
random tiling and some of its sub-ensembles to model decagonal quasi-crystal structures 
resulting from a molecular dynamics simulation. 

One of the main goals in statistical mechanics is the calculation of the partition sum 
which, for random tilings, is the weighted sum over all possible tiling configurations. 
To be able to enumerate all possible tilings we use the transfer matrix method, for 
which it is convenient to transform the tiling to a model on the square lattice. This 
transformation is depicted in figure ^ The triangles always come in pairs which are 
denoted by ti, . . . , ts corresponding to their five different orientations. Similarly the five 
different orientations of the rectangles are labeled by r^. The decorations on the deformed 
tiles are such that continuity of the decorating lines is equivalent to the restriction that 
the tiles fit together without gaps or overlaps. In this way every decorated tiling of 
the lattice with the deformed tiles corresponds to an allowed tiling of the plane by the 
original rectangles and triangles. Using this transformation every vertex of the tiling 
falls on a vertex of the square lattice, though the reverse is not true. The partition sum 
of the deformed model on the lattice is now defined by 

C i=l 

where the sum is over all possible configurations C. For a given configuration C, we 
denote the number of deformed rectangles r^ by iVrj(C) and the number of deformed 
triangles t, by Nt^{C). The partition sum is thus a weighted sum where each deformed 
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triangle tj and each deformed rectangle li is assigned a weight tj and rj respectively. The 
partition sum (|l]) is equal the partition sum for the tiling, provided that one chooses 
the weights of the deformed tiles properly |]1T|. The latter is a consequence of the 



transformation which changes the areas of the various tiles differently. 



tl 





Figure 1. Tiles of tlie the tenfold tiling with rectangles and triangles and their 
deformations to the regular square lattice. The five different orientations of the pairs 
of triangles and those of the rectangles are denoted by ti and r^ respectively. 



The definition of the transfer matrix T of the deformed model is the obvious one 
on the square lattice. The horizontal edges of the square lattice can be in one out of 
five possible states, represented by the presence or absence of the decorating lines. The 
transfer matrix elements Tij between two sequences i and j of such horizontal edges is 
equal to the product of the weights of the deformed tiles that fit in between i and j. 
As is well known, the partition sum per row in the thermodynamic limit is given by the 
largest eigenvalue of the transfer matrix. The free energy is then given by the logarithm 
of this largest eigenvalue. 

Recall that according to the random tiling hypotheses the entropy per area may be 
written as 

where the Ij are the quadratic phason strain invariants for the 3x2 phason strain 
tensor. The goal of our work is to calculate the residual entropy a^.^ and the phason 
strain elastic constants Kj. In section ^ we derive the Bethe Ansatz (BA) equations that 
diagonalise the transfer matrix for this model and which already give a huge reduction 
of the numerical problem. It turns out that for this tiling the method of Kalugin ^ is 
applicable to solve these BA equations . This then enables us to calculate the residual 
entropy for this tihng exactly. The calculation of the elastic constants poses some 
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problems and in this paper we will concentrate on the maximum only. It will be shown 
that the maximum of the entropy per vertex of this decagonal random tiling model is 
given by 

= jflog5-2v^logrV (3) 



2. Degrees of freedom 

In principle there are fifteen partial densities for this tiling, corresponding to the ten 
different orientations of the triangles and the five different orientations of the rectangles. 
Since the triangles always occur in pairs, we are left with ten degrees of freedom. One of 
these is removed by the fact that the total area is constant. Furthermore, there are three 
nonlinear geometrical constraints, so that the phase space of this random tiling is six- 
dimensional. The geometrical constraints are derived in the following.^ The first thing 
to note is that the triangles can be viewed as domain walls between patches consisting 
solely of one type of rectangle. This phenomenon is similar to what happens in de 
dodecagonal square-triangle and the octagonal rectangle-triangle tilings. In our choice 
of decoration in figure |I] it is easily seen that between patches consisting only of the 
tile 12 there are three different types of domain wall, denoted by the solid, dotted and 
dashed lines. Two types of domain wall run from bottom right to top left. They are 
drawn as solid and dashed lines and we denote their number by rii and n2 respectively. 
The other type of domain wall runs from bottom left to top right and is drawn as a 
dotted line. Their number is given by m. We denote the average number of triangles 
and rectangles per layer by nt^ and respectively. It is then easily seen from figure |l] 
that 

ni = K'^ta + n-tj + rir^ + ^r^, 

n2 = lnt^+n,^+n,,, (4) 
m = |(nt2 + rit^) + n^^ + n,^. 

Let pi be the number of layers such that each dotted line crosses each dashed line once. 
The number of such crossings in this patch of size piN is then 

nam = pi{n,^ + n,^ + |ntj. (5) 

Similarly, if p2 is the number of layers such that each dotted line crosses each solid line 
once, and q the number of layers such that each dashed line crosses each solid line once 

f This argument makes use of the periodic boundary conditions imposed by placing the lattice model 
on a cylinder. We believe that for free boundary conditions modified versions of such constraints hold, 
involving the configuration of the boundary. 
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we find 

mm = p2 (n,3 + + intg ) , 



nin2 = q{n,^ + ^15)- 

The numbers pi, p2 and q can also be calculated in another way. For this, it is convenient 
to introduce the average shift per layer (si, S2 and Sm) of each domain wall. These are 
given by 

51 = -1 + —Kg - n^g + ^ritj, 

Til 

52 = -1 + — - ^1.4 - hrit^), (7) 

n2 

From the condition that each dotted line must cross each dashed line once in a patch of 
size piN it then follows that in pi layers the average relative shift of the domain walls 
must be equal to N. The same arguments applied to the other cases then gives 

N = Pl{Sm - S2) = P2{Sm - Si) = g(si - S2). (8) 

Putting all these equations together, using the fact that the system size N = ni + n2 + 
m + 2(nr2 ~ ''^rs) one finds three independent relations among the tile densities. These 
can be rewritten as follows 

+ 2(nt3 +2(nr3 +nrJ)K, +nrj, (9) 



nt, (nt4 + 2 (^1.3 +n,^))+ An,^ (n,,^ +n,^), (10) 

which are symmetric under the simultaneous exchange of indices 2^5 and 3 4, i.e. 
the mirror symmetry in the y-axis. Other relations may be obtained from and (p!0|) 
by applying a rotation over 27r/5, i.e. shifting each index by one. Only three of them, 
however, are independent. 

At the symmetric point each orientation occurs equally often, so we have n^- = -g^tri 
and rir^ = ^n^ect, where ntri and nrect are the average total numbers of triangles and 
rectangles per layer. It follows then from and the expression for the system size that 
at the symmetric point ntri = lOr'^^A^ and n^ect = ^t~^N. Using (H) we find that these 
values correspond with the following numbers of domain walls. 



rn2 = rii = m = T "^N. 



(11) 
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3. Bethe Ansatz 

In this section we derive the Bethe Ansatz equations for the lattice model to diagonalise 
the transfer matrix. Again we make use of the fact that the triangles can be viewed 
as domain walls. Since these domain walls persist through the lattice we can think of 
them as being trajectories of three types of particles, where each number of particles 
is conserved. The transfer matrix T, acting in the upwards direction, can be thought 
of as an evolution operator for these particles. Two types of the particles are then left 
movers whose trajectories are given by the solid and dashed lines. We call them of type 
1 and type 2 respectively. The other type is a right mover and its trajectories are drawn 
as dotted lines. Because the number of particles of each type is conserved, T is block 
diagonal in the particle numbers. In the sequel we shall diagonalise T in each block 
separately by using a nested coordinate Bethe Ansatz [|l^ . 

A state a on a row of the lattice can be specified by the positions yi, . . . ,ym of the 
right movers and by the positions Xi, . . . ,x„ of the left movers with the specification 
that the lines ^l, . . . , im at positions Xj^, . . . , are of type 1. Elements ip{a) of an 
eigenvector of T thus can be written explicitly as . . . , i-m . . . , x^] • • • ? Vm)- 
The state which has no particles, i.e. the one with only the rectangles T2, is called the 
pseudo- vacuum. We make the following Ansatz for the form of the eigenvector, 

n m ni ^ ' ' 



EE^(r)5(/^)n4:n<n 



pb 

vr,p a=l 6=1 c=l L r=l 



(12) 



where the sum runs over all permutations /i = (/ii, . . . , of the numbers 1, . . . , ni, all 
permutations vr of the numbers 1, . . . , n and all permutations p of the numbers 1, . . . , m. 
The coordinates Xa and i/b enter the Ansatz (|l2|) as powers of complex numbers, Zi and 
Wj respectively, that are to be determined later. Factors in the eigenvector due to the 
order of the different types of right movers are given by the expression between brackets 
in ([T^) . This is the so called nested part of the Ansatz which is a generalisation of the 
simple power for the coordinates. The index ic may be seen as the 'coordinate' of right 
movers of type 1 relative to all right movers. This coordinate ic determines the upper 
limit of a product over a complex valued function u which may still depend on the 
numbers 2j. The Ansatz (|1^) further contains an adjustable constant d whose presence 
will become clear below. 

The amplitudes A depend on the permutations tt and p and on the configuration 
of the left and right movers. These together are coded in a vector F in the following 
way. Let p be the vector of coordinates Xi and i/j of all domain walls, ordered so that 
Pj < Pj+i- The entries of F are the elements of the permutations vr and p. The order of 
succession in F of elements taken from vr and p matches that of the elements of x and 
y respectively in p. The amplitudes B only depend on the permutation p. 
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If all the domain walls are separated the action of the transfer matrix is just a shift 
of each line to the right or to the left. The eigenvalue of T corresponding to the vector 
(|1^) is therefore given by 



A 



^iN-n-m)/2,2n-2ni,2m 



(13) 



a=l 



b=l 



where N is the size of the lattice. At places where different domain walls are close 
together, the action of T is not given by a mere shift of all domain walls. Whenever 
there is a right mover just in front of a left mover and there are no other neighbouring 
walls, two things can happen. Either the right mover jumps over the left mover, which 
does not move, or the left mover jumps over the right mover, see figure 0. These 



Figure 2. Two-particle collision diagrams. Upper line: Two types of collision of walls 
of type 1 and 3. Second line: Two types of collision of walls of type 2 and 3. 



exceptions imply the following relations for ([T^) to be an eigenvector of T. 



and they put the following constraints on the weights 

r^d = t\ and r^t\t1 = ritgt^. (15) 

The reason for the constant d in the Ansatz ([T2|) is now clear: if it were omitted there 
would be an additional constraint between the weights and t^. Scattering processes 
involving three particles of only two different types give the following restrictions on the 
amplitudes A, 

A{..., Ilk, TTfc+i, . . .) _ ^^TTfc+i Pi, pi+i, . . .) _ Wp^^^ 



A(. . . ,7rfc+i,7rfc, . . .) z^^ ' A{. . . , pi+i, pi, . . .) Wp^ 

From these processes it also follows that the amplitudes B obey the equation 



(16) 



B{...,pp,pp+i,...)z^^{u{pp+i,7rk)-u{pp+i,7ik+i)) = 0. (17) 

By the sum in (^) we mean a sum over both permutations of the numbers Pp and 
Pp+i- A similar notation is used in the next two equations, which result from scattering 



processes with three particles which are all of a different type. 

A{...,7rk,7rk+i,...)z^^^u{iJ,p,'Kk) = 



Af \ ^^^^^^ -1 I ^1*^4 2 \ ns^ 



I \ I '2'3'5 -2 i-lM 2 /1n^ 



Equation (|TBp is satisfied (after A is eliminated using (|TB|)), if u is of the form 

with any complex number f^^. Substituting (|T9|) into (p!7| ) it follows that the amplitudes 
B fulfil the relation 

B{. ■ . ,iJ,p, Atp+i, ■ ■ ■) ^ _^ ^20-) 
. . , fip+i, fip, ■ ■ .) 

Using periodic boundary conditions, it follows from the form of the eigenvector (|T2|) and 
the relations (p!^), (|1^) and ( |20D that the complex numbers Zi, wj and Vk should obey 
the following Bethe Ansatz equations 



■)'"-' fi (^) n ih^^', + ^^i'y^A . (21) 



k=i ^^i' j=i V2 ^1 / 



1=1 



i2 \ ^ n / ^ ^ ^ +2, 



To rewrite the BA equations in a more suitable form we introduce the following variables, 
( r4\t\ 2 , / ^i^4t2 \ ^'^ -2 



\r1r2r3r5y 

1/2 



-1 / \ 



The BA equations then become 



m 



(_^^.)(7V+n+™)/2 ^ (_)-lCne-'/' m'^^f'"Wl^ - ^.), (25) 

i=l k=l i=l 

n m rti 

gA.+n+.)/2 ^ n(-^.)-^/^ n%-^ X 

k=l j=l 1=1 

m ni 

- ^j) Uiui - ^,){Ui + C'), (26) 

j=l 1=1 
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u- = i-r-'El[iui-mui + ^r'), (27) 



i=l 

where C, D and E are given by 

N/4: 

I I -I I A I, ^ \ 

c 




\r2r3r5J 



N / N , 



The eigenvalue in terms of these new variables is then given by 

In summary we have shown in this section that the Bethe Ansatz equations (]25|), 
(p6|) and ( pTj) with the definitions (p8D diagonalise the transfer matrix T for arbitrary 
choice of the weights except for the one constraint rstftg = rit|t|. It thus follows that 
the point of maximum symmetry, where rj = r for i 7^ 2 ^ and tj = t, is included in 
the spectrum obtained by this Ansatz. The eigenvalue of T in terms of solutions of the 
Bethe Ansatz equations is given by (p^). 

4. Integral equations 

To calculate the entropy, we put tj = 1, r2 = e'^^ and ri = = = = e^^ . The 
chemical potentials /zi and /i2 have to be adjusted such that all configurations in the 
original undeformed tiling model are weighted properly. The difference between /ii and 
/i2 compensates for the fact that the area of the transformed tile r2 is twice that of the 
other transformed rectangles. If we want to weight the original rectangles in the random 
tiling equally it follows that they must satisfy [|r^ 

/il - /i2 = N'^ log Amax- (30) 

One finds numerically that the maximum of the entropy of the tiling model is at the 
point of maximum symmetry, i.e. where all configurations are weighted equally, in 
agreement with the first random tiling hypothesis [14|. According to ( [TTl) this is in the 



sector rii = m = T "^N, n2 = r '^A^, corresponding to tile fractions nrect = |t and 
= lOr-^A^. 

It is also observed numerically that each of the roots ui approximates one of the 
roots C,i in exponentially good precision, see figure ^. For reasons that will become clear 

f We need to tune r2 to compensate for the change in area induced by the transformation to the lattice. 
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Figure 3. Distribution of roots for the largest eigenvalue (N — 89, ni = m = 34, = 
21). The left curve corresponds to the roots f/'i and the outer right curve corresponds 
to some subset {^i} of the roots ^i. The inner right curve actually consists of two 
curves corresponding to the roots and the subset of the roots complementary to 
{^i}. Note that the inset is not on scale. 



below, we define the roots now as the subset of the roots not approximated by one 
of the Ui and we introduce the following notation 

n2 ra n\ 

H = V^i\'\ H = V^^-^^f'\ Su = ViA'"- (31) 
i=i j=i i=\ 

By writing ^ = + for those roots ^ that are approximated by one of the ui and using 
the above abbreviations, equation (^) splits up in the following two sets of equations, 

ra n\ 



i=i 

711 

^{ui -Uk- £k) (1 + uT^iuk + Ek)'^) ■ (33) 
1=1 

Similarly, substituting = u + e in (|27|) results in 

i-r~'E''ur = f[iui-Uk- Sk) (1 + ui\uk + Sk)-'] 

k=l 

n(«^-^«)(^/+c')- (34) 



X 



i=l 
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In the thermodynamic hmit all Ek vanish exponentially in A^. Equation p^ ) can then be 
used to remove the product over the variables u in (|33D and we arrive at the following 
equation which approximates the original BA equation (|33|), 



u 



{N+ni-n2+m)/2 



X 



-N\ 



(35) 



i=l 



Consider now the BA equation (|25D as a function of ipj. Taking the logarithm on both 
sides of ( p5D we define the function by 



F4z) = log(- 



N + rii + n2 + m 



■ n2 ni 
.i=l 1=1 



- jNni + |(ni + n2)/i2 - log 



(36) 



so that ^F^{ipj) = 0. Similar functions F^{z) and ^^(2;) are defined by taking the 
logarithm of equations ( p2D and ( |35D respectively. The BA equations (p5D, (|32| ) and (p5|) 
are then equivalent to 



i(A^ + rii + ria + m)F^{%l)k) = 27ri4, 
i(Ar + 3ni + 712 + m)Fg(a.) = 27riJfc, 
i(A^ + rii - 722 + m)F„(Mfc) = 27rii^'fc, 



(37) 



where Ik, Jk and Kk are either integers or half- integers. From numerical calculations it 
follows that the numbers 7^, Jk and Kk for the solutions of the BA equations for the 
largest eigenvalue are consecutive, more precisely 



(38) 



h 


= |(m + 1 


-2k), 


{k = 


1,.. 


. ,m) 


Jk 


= 1(^2 + 1 


-2k), 


{k = 


1,.. 


• ,^2 


Kk 


= U^l + 1 


-2k), 


{k = 


1,.. 


■,ni 



We will assume that ( p8| ) holds in the thermodynamic limit. It is for this reason that 
we introduced the variables ^ instead of ^. The function analogous to (|36| ) defined from 
(^) does not take consecutive multiples of 27ri when evaluated at the roots ^k- It is 
clear from (|38D that the derivatives / of the functions F are up to a factor precisely 
the densities of the Bethe Ansatz roots. They allow us to transform the sums in the 
logarithmic form of the BA equations into integrals. For brevity we define 



a 



1 + 3Qi + Q2 + Qn 
1 + Q 



P 



1 + Qi - Q2 + 

1 + g 



(39) 
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where Q = Q1 + Q2 + Qm- Taking into account the root distribution ( pHD we thus arrive 
at the following integral equations for the functions /, 



z 2711 J^i z — ^ 27ri Jui 



z — u 



(40) 



z 2ma 
1 

+ 



'"^ Mild. + !^:!Mz!^d. 



Z — M 



Z — M 



2nia JiPi z — ip 



d^, 



(41) 



fu{z) 



a 



z 2^13 
1 

+ 



hii) 



-1^ 



-de 



d^, 



(42) 



2mf3 Ji>i z — ip 

where the integrals are taken along the locus of the roots C,, 'ip and u and —ip'^ 
and —u^^. These integration contours in the complex plane are schematically shown in 
figure |. 




Figure 4. Schematic picture of the solution curves of the BA equations. 



5. Monodromy properties 

From the integral equations it is immediately seen that the integration contours 
are cuts in the complex plane of the functions f^{z), f^{z), fu{z), z^"^ f^{—z~^), 



13 



z'"^ f^{—z~^) and z""^ fu{—z~^). Each of these functions has jumps across some of the 
cuts whose magnitudes are given by hnear combinations of the above six functions. 
The analytic continuation across the cuts of these functions is now determined by 
compensating for the jump and can be written in terms of monodromy operators, 
one for each cut. The analytic continuations across the different cuts of a function 
G[z) = Yji=s^,-^,u{.'^ifi{.z) + fi{—z~^)) are given by the following matrices which act 
on the vector a = (a^, a^, Qu, b^, b^, bu), 

=1- a-'E2i - p-'E23, 
r^-i =1- a'^E54 - r^E^e, 

=I + aEi2-af3-\Ei3 + Eie), 
r^-i = J + aE45 - a(3-\E^3 + E^e), 

=1 + (3E32 + f3a~\E3i + Eu), 
r„-i = / + pEe, + Pa-\Ee, + E^^), 

where / is the 6x6 unit matrix and Eij is the matrix with a 1 at the entry ij and O's 
everywhere else. As a result, several linear combinations of the six functions, following 
from these monodromy operations, correspond to different sheets of the Riemann surface 
of one function. 

The continuation of a function G{z) starting at and returning to the origin across all 
of the six curves is given by the monodromy operator F = F^-ir^-ir^r^r^r^-i which 
has the property that = 1. This means that if the curves in figure ^ have the same 
endpoints, say b for the common endpoint in the upper half plane and b* for that in the 
lower half plane, that function is single valued in terms of the variable s: 

= : — T , z(s) = b — (44) 

^ ^ \l-zb*-^J ^ ' l + bb*-^s^ ^ ' 

When the BA-curves have the same endpoints, F is the only non trivial monodromy 
operator and the Riemann surface of G{z) breaks up into infinitely many disconnected 



parts each with only five sheets. These are all mapped onto the plane by (|4^), from 
which it is easily seen that, apart from the branch points z = b and z = b*, each point in 
the z-plane has five images on the s-plane. As can be seen from figure |^, for the largest 
eigenvalue the curves do indeed have the same endpoints, in particular they meet at 
b = —b* = i. The situation is rather similar to what happens in the square triangle 
0] and the octagonal rectangle triangle tiling 0, where also the BA curves close, but 
there the monodromy is of order 6 and 4 respectively. Although we do not understand 
the deeper reason for this, it will become obvious in the sequel that the order of the 
monodromy must relate to the symmetry of the tiling to produce the right quadratic 
irrationalities. 
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In the following we show that the closing of the curves enables us to calculate the 
largest eigenvalue explicitly. Note that the six curves can only have the same endpoints 
if 6 = —b^^ = i, see figure ^, in contrast to the square triangle and the octagonal 
rectangle triangle tiling, where the common endpoints need not lie on the imaginary 
axis. We will, however, use the notation b = i|6|e^''' for the common endpoint of the 
curves, for it may turn out in the future that for infinitesimal values of 7 still something 
can be said. 

Aside from the cuts, the form Gdz has only simple poles ai z = and at 2; = 00. 
The poles and their residues of the single valued function are thus given by the poles 
and residues on each of the sheets. The functional form of G{z) on other parts of the 
sheet and on the other sheets is directly determined by the monodromy operators (|43|) . 
In particular, the value of the residues can then be read off from the integral equations 
(1^), (^Tj) and (|42|) . We choose G{z) to be equal to f^p^z) near z = on the sheet 
with Si = e''^/^. Its form near z = on the sheets with S2k-i = e*^^'^"^-*''^/^ is obtained 
by applying T^. The functional forms of the function G{z) near z = cxd on the sheets 
with S2k = e^'*^^'^"'''^/^ are obtained by analytic continuation, i.e. by applying r^r^r^-i 
on each of the functions corresponding to S2k-i- The poles s„ and residues Rn thus 
obtained of this function G are shown in table |l|. 



Table 1. Poles and residues of the form G{z)dz. The two columns on the left list 
the poles of G{z)dz on the z- and on the s-plane respectively. The third column gives 
the functional form of G{z) near those poles and the fourth column the value of the 
residues of these poles. 



z 


Sn 


G 


Rn = Res- 


i(Gdz) 





g7ri/5 




1 




00 


g2i(7r-7)/5 


U{z)+PUz) 


-2(1+Q2) 
1+Q 







_g-27ri/5 


f3Uz)-az-^f^{-z-^) 


2 

1+Q 




00 


_g-i(7r+27)/5 


-az-^f^{~z-^) + z-^U{~z-^) 


-2Qi 
1+Q 







-1 


z-'Ui-z-') 


1+Q™-Qi 
1+Q 


-Q2 


00 


_gi(7r-27)/5 


z-^Ui-z-^) 


1 







_g2i7r/5 


z-^U{-z-^)+(3z-^fu{-z-^) 


-2(1+Q2) 
1+Q 




00 


g-2i(7r+7)/5 


(3z'^U-z-^)~aJ^{z) 


2 

1+Q 







g-iir/5 


-af^{z) + f^{z) 


-2Qi 
1+Q 




00 


e-2i7/5 




i+Q™-Qi 
1+Q 


-Q2 
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The form Gdz is now uniquely determined by its poles and residues, 

10 p 

Gdz = g -^ds. 



n=l 



S — Sr 



(45) 



Recall that we are working with approximated BA equations which are valid for the 
largest eigenvalue, for which Qi = Qm = Q2 = and thus Q = 1. We do not 
know at the moment if the approximation is valid for other sectors as well. The closing 
of the three curves at 6 = i (7 = 0), however, is a fact that only holds for the largest 
eigenvalue. The above arguments for calculating Gdz therefore are valid for this sector 
only, for which Gdz is given by 



Gd^=(se-2-/^ + s-V'-/^)r-i— . 



(46) 



6. Solution of definite integrals 



In this section we calculate some definite integrals of the form (ffSp, resulting in an exact 
expression for the maximum of the entropy. To remove singularities in the expressions 
to follow, we introduce the following forms 

R2n-1 



Qndz 



-ds. 



(47) 



S — S2n-l 

In the 2;-plane we can calculate the following integrals using the fact that 6 = i is a 
solution of the BA equations and that we have a precise control of the singularity at 



2; = by using the mapping (|4^) to the s-plane 



J' 



2,0 



J- 



3,0 



J. 



4,0 



J' 



5,0 



iUiz) -gi)dz 



3? hm 

g_,giTr/5 



F^{z{s)) - log(s - e 



i7^/5^ 



log(5/2) - N^^ log A + /ii- - fi2r 



-2 



3? 



{f3fu{z) - az~'^f^{-z'^) - g2)dz 



log(5/2) - 2N-^ log + fxi-r-^ - /i2r" 



3? 



3? 



3? 



iz~'^Ui-z-^) - 93)dz 



log(5/2) + log A - log - /i2r-^ 

\z-'Ui-Z-') + I3z-'f^(-z-') - g,)dz 

2r-i log(5/2) + /iiT-i - fi22r-\ 



'-af^{z) + U{z)) -g5)dz 



1 1 
log(5/2) + fii^r-^ - f^2r'^ = ^ ^o- 



(48) 

(49) 

(50) 
(51) 

(52) 
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The same integrals on the s-plane restricted to the symmetric point are, using (^), 

10 

Ji,o = ^nlog h"'''^ - Sn\ = r logr - logsin(27r/5), 

n=2 

^2,0 = rlogr - logsin(27r/5), 
-^3,0 = -r"Mogr + r~Mogsin(27r/5), 
J4,o = -21ogr + 2t~^ logsin(27r/5), 
J5,o = -T^^ logr + r"^ logsin(27r/5). 

Equating both sets of integrals gives the following solution 



(53) 



4r~^s^ = /ii - (1 + r-^)/i2 = A^~^ log A - r^^/xa 

1/5^ \ 
= 2 -2v^logrj . 



(54) 



Note that this solution precisely corresponds to the maximum of the entropy for the 
original tiling model according to (|30|) . The total number of rectangles per site, 
Qrect = n^-ect/N, ou the symmetric point is |t~^, so that the entropy per site is given by 

aj^ = log A - |r-^(/ii + 4/i2) = It-\N-' log A - r-V2). (55) 

The number of vertices per site is given by 

nv = iQtri + Qrect = ^'^2 + T'') = ^T'^ (56) 

Thus, the entropy per vertex is finally given by 

a, = log A - r-V2 = ^ (log ^ - 2^5 logr j . (57) 



7. Conclusion 



In this paper we showed that a decagonal random tiling model of rectangles and 
triangles is solvable using the Bethe Ansatz technique. We derived the Bethe Ansatz 
equations that diagonalise the transfer matrix for this model. These equations contain 
all information about the model and they in principle present a huge reduction of 
computational problems concerning the system size. For this tiling model however, some 
of the roots of the BA equations almost coincide, which makes it difficult to extract 
high precision data. On the other hand, it enabled us to write down approximate 
BA equations which are exact at the symmetric point in the thermodynamic limit. 
Using these equations we were able to find an exact expression for the maximum of 
the entropy. The validity of the approximation outside the symmetric point still has to 
be further investigated. We hope to be able to find analytic expressions for the phason 
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elastic constants as well, although in contrast to the solutions of the dodecagonal square- 
triangle 0, and the octagonal rectangle triangle 0] it is not apparent if our solution 
can be extended off the symmetric point. Recently, Kalugin |T^ showed for the square 
triangle model that also critical exponents may be calculated exactly from the BA 
equations. 
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